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In this paper we categorize dynamical regimes demonstrated by star-like networks with chaotic 
nodes. This analysis is done in view of further studying of chaotic scale-free networks, since a 
star-like structure is the main motif of them. We analyze star-like networks of Henon maps. They 
are found to demonstrate a huge diversity of regimes. Varying the coupling strength we reveal 
chaos, quasiperiodicity, and periodicity. The nodes can be both fully- and phase-synchronized. The 
hub node can be either synchronized with the subordinate nodes or oscillate separately from fully 
synchronized subordinates. There is a range of wild multistability where the zoo of regimes is 
the most various. One can hardly predict here even a qualitative nature of the expected solution, 
since each perturbation of the coupling strength or initial conditions results in a new character of 
dynamics. 

PACS numbers: 05.45.-a, 05.45.Xt, 05.45.Pq, 89.75.Hc 


I. INTRODUCTION 

Complex dynamical networks with scale-free coupling 
structure attract a lot of interest as models for a large va¬ 
riety of natural systems [1, 2]. The name “scale-free” for 
these networks appears because their node degree dis¬ 
tributions have power law shapes. As a result a small 
number of nodes hold a major balk of links while the 
rest of nodes have few connections [3] . 

One of the main questions arising in studying of these 
networks is the type and conditions of synchroniza¬ 
tion [1, 4-6]. The synchronization can be full or only 
the phases of node oscillators can be synchronized; it 
can involve the whole bunch of nodes or the nodes can 
form synchronized clusters [2, 7, 8]. Papers [4, 5, 9, 10] 
are devoted to an analysis of synchronization conditions, 
works ]11, 12] investigate the formation of synchroniza¬ 
tion clusters while in Ref. ]13] the impact of presence 
of a leader on the cluster synchronization is recovered. 
Paper ]14] investigates so called remote synchronization 
when nodes can get synchronized even being connected 
indirectly though intermediate ones. Also this regime is 
studied in Refs. ]11, 12] being called driven synchroniza¬ 
tion. 

Authors of Ref. ]15] study scale-free networks with 
fractional order oscillators. Paper ]16] investigates the 
control of a scale-free dynamical network by applying lo¬ 
cal feedback injections to a fraction of network nodes. 
Covariant Lyapunov vectors ]17] and their nonwandering 
predictable localization is studied in Ref. ]18]. 

The feature specific for scale-free networks as well 
as for other complex dynamical networks is multistabil¬ 
ity ]19-21]. It is well known that the dynamics of multi¬ 
stable systems can be amazingly rich [22, 23]. It occurs 
when the number of attractors is very high, and their 
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basins have fractal boundaries that are highly interwo¬ 
ven. In this case the dynamics is very sensible to the 
initial state: even tiny perturbation results in arriving at 
new regime. Moreover, the ranges of existence of par¬ 
ticular attractors can be narrow so that the qualitative 
behavior of the system can change dramatically when its 
parameters are slightly varied ]22]. For scale-free net¬ 
works this type of behavior was reported in Ref. ]18]. 
We suggest to refer to this type of dynamics as wild mul¬ 
tistability to distinguish it form the plain case, when one 
can easily locate basins of the required attractors and put 
the system there to observe the expected behavior. 

Wild multistability as well as other dynamical phenom¬ 
ena demonstrated by dynamical networks still requires 
an exhaustive study. In particular the interesting ques¬ 
tion is to reveal what features are the same with other 
systems and what additionally emerge due to the strong 
inhomogeneity of networks. For scale-free networks the 
this study can be started from the consideration of their 
simplest and regular representatives which are star-like 
structures. These structures consist of the hub node con¬ 
nected with all other nodes and subordinate nodes that 
have only one connection with the hub. Star-like struc¬ 
tures are the main motifs of scale-free networks, their 
building blocks. 

Chaotic synchronization of oscillator networks with 
star-like couplings is considered in Ref. [24]. Formation 
of synchronized clusters in such networks is studied in 
Ref. ]25] and a sufficient condition about the existence 
and asymptotic stability of a cluster synchronization in¬ 
variant manifold is derived. Paper ]26] considers phase 
synchronization of the subordinates when the hub is not 
synchronized with them. This is called remote synchro¬ 
nization. 

Thus, to pave the way to understanding the dynamics 
of scale-free networks, one has to reveal first the details 
of dynamics of systems with star-like coupling. This is 
the main motivation for the present paper. We consider 
Henon maps. This map is a canonical model that is suf- 
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ficiently simple on the one hand and exhibit the same 
essential properties as much more “serious” chaotic sys¬ 
tem on the other hand [27]. In particular, this map is 
time-reversible, and the coupling between network nodes 
is introduced in a way preserving this property. Alto¬ 
gether, we build a catalog of regimes of time-reversible 
star-like network of Henon maps observed at different 
coupling strengths. The most of the regimes are found to 
be independent on the network size. However there is an 
area of wild multistability where the dynamics does de¬ 
pend on the number of nodes. Many different attractors 
exist here within a narrow ranges. As a result, one can 
hardly predict even a qualitative nature of the expected 
solution, since each perturbation of coupling strength or 
initial conditions leads to a new character of dynamics. 

The outline of the paper is the following. First we in¬ 
troduce the model system. Then we discuss the numeri¬ 
cal criteria for detecting various synchronization regimes. 
After that the registered regimes are represented and dis¬ 
cussed. Finally, the ranges of stability of some of the 
regimes are derived. 



Figure 1. (color online) The star-like network, A = 4. 

where ^{t) G (—10“^^, 10“^^). Notice that the noise is 
added only to fix the drawback of the numerical model. 
Since its amplitude is not much higher then the machine 
epsilon for employed variables, « 10“^®, it is indistin¬ 
guishable in the resulting time series. Thus, we empha¬ 
size that the study of noise influence on the network lies 
outside the scope of our paper. 


II. MODEL SYSTEM 


III. CRITERIA OF SYNCHRONIZATION 


We consider a network of Henon maps introduced in 
Ref. [18] as a generalization of the Henon chain from 
Ref. [28] : 


Xn{t -k 1) = a - [Xni,t) + ehnit)]'^ + yn{t), 

yn{t + 1) = j3Xn{t), 

N N 

hn{t) ='^-^Xj{t) - Xn{t), kn = '^ajn- 

.--I 


( 1 ) 


Here N is the number of network nodes, f = 0,1, 2 ... is 
discrete time, a„j G {0,1}, a„j = ajn, = 0 are the 
elements of the N x N adjacency matrix A, and is 
degree of the nth node, i.e., the number of its connections. 
a = 1.4 and /? = 0.3 are the parameters controlling local 
dynamics, and e G [0,1] is the coupling strength. Recall 
that the Henon map is time-reversible. The coupling is 
introduced in a way that preserves this property. 

We consider star-like networks, see Fig. 1. Let e,n,Tn 
be a unit mxm matrix, 2^,1 be a m x n matrix of zeros, 
and dm,n be a to x n matrix of ones. Using this notation 
the adjacency matrix can written in block form as 

di,M\ (2) 

\dM,i zm,m/ 

where M = N — 1. 

In numerical simulations some regimes was found to 
emerge outside the theoretical ranges of stability, com¬ 
puted in Sec. V. In particular this is the case for full 
chaotic synchronization. This spurious stability occurs 
due to round-off errors in numerical model. To eliminate 
it a very weak noise is added to system variables: 


Xn(t) —>■ Xn(t) -k yn(t) —> yn{t) + + 1), (3) 


Coincidence of local extrema of two discrete time sig¬ 
nals can be treated as their phase synchronization while 
the relative frequency of the coincidences quantifies the 
degree of the synchronization [11, 12[. Below we consider 
a generalization of this approach taking into account both 
dip-dip and peak-dip coincidences. 

Let us consider the nodes m and n, m < n. Given a 
starting time to and the interval T, count at to < t < 
to +T the numbers Vm and of local minima of Xm (t) 
and Xn{t), respectively, and the number i/^n of simul¬ 
taneous minima of Xm and Xn- In addition we count 
the number of times jimn when the local minima of Xm 
occur simultaneously with the local maxima of Xn, and 
the number of times ynm when the minima of coin¬ 
cide with the maxima of Xm- Note that since m < n, 
fimn can be treated as elements of the upper triangle of 
a square matrix and form its lower triangle. Then 
the dip-dip distance can be computed as 


d 


d 

m-n 


1 - Vmn! max{z/m, i/„} if max{r'm, Vn} > 0, 

0 if max{r'm,^'n} = 0. 

( 4 ) 


The second choice covers the situation when both x^ 
and Xn do not oscillate at all. This value is introduced 
in Refs. [11, 12] as phase distance. Analogously we can 
define peak-dip distance as 


dP = 


1 - l^mn/max{Vm, t'n} if max{Vrn, l^n} > 0, 

0 if max{z/m,r'„} = 0. 

( 5 ) 

Two nodes are did-dip or peak-dip synchronized on the 
interval T if = 0 or = 0, respectively. 
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The full synchronization can be detected nsing 
dL-n = Var{a;m(t) - Xn{t) \to <t <to+T}, (6) 

where Var stands for variance. Moreover we are going to 

detect complementary synchronization via the vanish of 
■ 

^m-w 

d^-n = Var{a;^(t) + a;„(t) | to < * < + T}. (7) 

The synchronization is actually registered when or 
are below the threshold 10“^^, which the level of 
the added noise, see Eq. (3). 

All of the above criteria depend on T. It determines 
the resolution of the detection procednre. If T —t oo the 
procednre responds only to fnll-time regimes. So it is 
preferable for T to be as short as possible. In this case 
in addition to full-time synchronization, one can detect 
synchronization windows as series of subseqnent inter¬ 
vals where d^_„, or vanish. The lower 

boundary for T is the average interval between peaks and 
dips: T has to be large enough to catch at least two-three 
of them. In simulations below we set T = 16. 


IV. CLASSIFICATION OF DYNAMICAL 
REGIMES 

A. Method of analysis 

Since it is unclear a priori what types of behavionr 
can be observed, the most reliable way is the visnal in¬ 
spection of time series corroborated by some appropriate 
characteristic nnmbers. 

First of all we are going to employ the first Lyapnnov 
exponent whose positive value indicates chaos, the neg¬ 
ative sign reveals periodicity and the zero means qnasi- 
periodicity^. One has to take many random initial condi¬ 
tions and compute Ai for each corresponding trajectory. 
The resnlting values are grouped very well near a few 
points representing different regimes. This approach is 
nsnally nsed for analysis of mnltistability [23] . 

Varions types of synchronons oscillations will be de¬ 
tected nsing criteria introduced in Sec. III. We take 
T — IQ and test at each snbsequent interval T if one 
of the fonr characteristic valnes vanishes for each pair of 
oscillators. If yes, we are inside a window of synchroniza¬ 
tion of a corresponding type. Then we find the largest 
lengths of the windows that are registered along the ob¬ 
servation time, which is tmax = 10®T. The full-time syn¬ 
chronization is registered if the corresponding window 
lasts during the whole observation time tmax- 


Table I. Stability ranges of e for regimes in Figs. 2 and 3 
obtained via step by step analysis of dynamics with Ae = 
0.01. More accurate values for FSC, OSD and FCP are found 
in Sec. V, see Eqs. (14), (21), and (26), respectively. 


NSC 

[0,0.11], [0.25,0.34] 

FPQ 

[0.12,0.14] 

DPP 

0.14 

FCP 

[0.15,0.24] 

FaPP 

0.15, 0.16 

FAP 

0.18, 0.82, 0.88 

FaPC 

[0.19,0.24] 

FSC 

[0.35,0.82] 

OSD 

[0.76,0.85] 

FNQ 

[0.86,0.88] 

FNC 

[0.89,1] 




B. The star with N = 4 

Consider dynamics of the network (1), (2) when the 
coupling e varies from 0 to 1. We have inspected all 
values of e S [0,1] with the step 0.01 for the smallest 
nontrivial star N = 4. Moreover when a regime of some 
sort appeared only on a single step, we also checked if it 
existed at least in a small vicinity of the corresponding e 
to make snre that this is a typical situation. To detect a 
mnltistability, for each e we tested at least 200 random 
initial conditions. The results are gathered in Figs. 2 and 
3 that are discnssed below. 

To facilitate the referencing we will label the regimes 
with three letters. The first two ones are an abbrevia¬ 
tion describing a type of synchronization between net¬ 
work nodes. In some cases we will snpply them with a 
snbscript to show a nnmber of involved nodes. The third 
letter indicates the character of oscillations: periodic “P”, 
qnasiperiodic “Q”, chaotic “C”. If there are several regimes 
with identical abbreviations we also add an index nnm¬ 
ber. The exception to this scheme is the oscillation death 
that will be labeled merely as OSD. 

The ordering of the regimes is shown in Fig. 3(a), and 
the boundaries of the regimes are collected in Tab. I. 
Fig. 3(b) shows Ai vs. e computed for 200 random initial 
conditions. 

• e = 0,0.01,..., 0.11: NSCl^. There is only one 
chaotic attractor here. This is confirmed by the coinci¬ 
dence in Fig. 3(b) of Lyapunov exponents, and by visnal 
inspection of time series whose typical example is shown 
in Fig. 2(a). The power law divergence of peak-dip syn¬ 
chronization windows between the hub and snbordinates 
and dip-dip windows for the snbordinates, see Figs. 4(a) 
and (b), indicate the intermittency appearing near the 
right bonndary of the considered range. 

• e = 0.12,0.13: FPQ^, Fig. 2(b). The quasiperiodic¬ 
ity is confirmed by the vanish of Ai in Fig. 3(b). 


^ In principle, a negative Lyapunov exponent can also be observed 
for so called strange non-chaotic attractors [29]. However, we did 
not encountered them in our case. 


2 

3 


No Synchronization, Chaos 
Full synchronization of the 
Peak-dip synchronization of 
Quasiperiodicity 


subordinates, = 0, 

’ sub-sub ’ 

the hub with them, dF u u 

’ nub-sub 


and 
= 0 , 
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-1-1-1-1-1-1-1“ 

(b) FPQ, e=0.12 


-1-1-1-1-1-1-1“ 

(c) DPP, e=0.14 


-1-1-1-1-1-1-1“ 

(d) FCP, e=0.15 



■ ■ ■ 




(e) F 2 PPi !£=015 


“I-1-1“ 


(f) F 2 PP 2 , e=0.16 


(g) FAPl,e=0.18 


(h) F 2 PC, e=0.19 




■ ■ ■ 




-1-1-1-1-1-1“ 

(i) NSC2, e=0.26 


“T-1-1-1-1-1-1“ 


“I-1-1-1-1-1-1“ 


-1-1-1-1-1-1-1“ 

(1) FAP2, £=0.82 


(j) FSC, £=0.5 


(k) OSD, £=0.81 




I I I 



(m) FNQ,'£=0.87 


“I-1-1“ 


(n) FAP3, £=0.^ 


(o) FNC, £=0.89 


(p) FNC (rich), £=0.97 



0 4 


12 16 20 24 28 t 


Figure 2. (color online) Various regimes of the star-like network (1), (2) with V = 4. See details in the text. 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 e 

Figure 3. (color online) (a) An ordering of the regimes shown 
in Fig. 2. See also Tab. 1. (b) The first Lyapunov expo¬ 
nents computed for 200 random initial conditions at each 
e = 0,0.01,0.02,...,! 


• e = 0.14: FPQ and DPP^, Fig. 2(c). Actually in 
DPP two subordinates are fully synchronized while the 


^ Dip-dip synchronization of the subordinates, = 0, and 


third one, the node number 4, almost coincides with 
them. The Lyapunov exponent corresponding to DPP is 
close to zero, Ai = —0.00030, and one can not distinguish 
it in Fig. 3(b). However, the inspection of data reveals 
that it is strictly smaller then the Lyapunov exponent 
for FPQ, so that despite of FPQ this regime is periodic. 
Moreover, we computed Fourier spectra (not shown), 
that again confirmed the quasiperiodicity of FPQ and 
the periodicity of DPP. 

• e = 0.15: FCP^ and F2PP1^. See Figs. 2(d) and 
(e), respectively. The dominating regime here is FCP. It 
exists within a wide range of e, see Tab 1. In Fig. 3(b) 
it corresponds to the lowest branch of the Lyapunov ex¬ 
ponents. The F 2 PPI has smaller basing of attraction 
and exists within a narrow range of e around 0.15. At 
e = 0.15 the Lyapunov exponents for these two regimes 
though different, are close to each other and barely dis¬ 
tinguishable in Fig. 3(b): for FCP Ai = —0.0487 and for 
F 2 PPI Ai = -0.0389. 

• e = 0.16: FCP and F 2 PP 2 . Here we have the sec¬ 
ond version of the periodic regime with full synchroniza¬ 
tion of two subordinates, see Fig. 2(f). Though it looks 


Peak-dip synchronization of the hub with them, = 0, 

Periodicity 

^ Full synchronization of the subordinates, = 0, and Com¬ 
plementary synchronization of the hub with them, = 0 , 

Periodicity 

® Full synchronization of two of the subordinates, ^ 2.3 = O 5 Peak- 
dip synchronization of the hub with this two, d ^_2 3 = Ihe 
third subordinate oscillates separately, Periodicity 
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like F 2 PPI, the closer inspection reveals that the period 
of F 2 PPI is 6 and the period of F 2 PP 2 is 20. The first 
Lyapunov exponent for F 2 PP 2 is Ai = —0.00191. 

• e = 0.17: FCP. 

• e = 0.18: FCP and FAP1~, Fig. 2(g). The Lya¬ 
punov exponents are Ai = —0.243, and Ai = —0.0360, 
respectively. 

• e = 0.19,..., 0.24: FCP and F 2 PC®, Fig. 2(h). 
The emergence of the second regime is clearly seen in 
Fig. 3(b), where the upper branch of positive Lyapunov 
exponents appears near 0.2. Note that Ai for F 2 PC is 
rather independent on e. 

• e = 0.25,..., 0.34: NSC2, Fig. 2(i). In the sec¬ 
ond version of these regime we again observe an inter- 
mittency. When e approaches the right boundary of 
the range the full synchronization windows demonstrate 
power law divergence, see Fig. 4(c,d). 

• e = 0.35,..., 0.75: FSC®, Fig. 2(j). In Fig. 3(b) this 
regime is represented by a perfect horizontal line. 

• e = 0.76,..., 0.81: FSC and OSB^o, Fig. 2(k). OSD 
has very small basin of attraction here. It becomes more 
or less easily detectable only at e = 0.81. There are two 
forms of OSD. The second one, not shown, is an inter¬ 
change of the values for the hub and subordinate node 
variables. Two forms have different negative Lyapunov 
exponents, so that there are three branches in Fig. 3(b), 
one positive for FSC and two negative for OSD. 

• e = 0.82: FSC, OSD, and FAP2, Fig. 2(1). Ac¬ 
cordingly, there are four values of Lyapunov exponents, 
see Fig. 3(b): Ai = —0.243 and —0.158 for OSD, Ai = 
-0.0492 for FAP2, and Ai = 0.419 for FSC. 

• e = 0.83: OSD. The Lyapunov exponents corre¬ 
sponding to the two forms of OSD becomes very close 
to each other, Ai = —0.192 and Ai = —0.182. Right 
after this point they merge, see Fig. 3(b). 

• e — 0.84, 0.85: OSD. Now both forms of this regime 
have identical Lyapunov exponents, see Fig. 3(b). 

• e = 0.86,87: FNQ'^^. This regime appears when two 
fixed points corresponding to OSD become unstable, see 
Fig. 2(m). The corresponding Lyapunov exponent is zero 
here, see Fig. 3(b). 

• e = 0.88: FNQ, FAP3, Fig. 2(n). The Lyapunov 
exponents are zero for FNQ and Ai = —0.0888 for FAP3, 
see Fig. 3(b). 

• e = 0.89,..., 1: FNC^^. Starting from this point the 
oscillations born from OSD fixed points become chaotic. 
Initially the oscillations have sufficiently small amplitude, 
see Fig. 2(o), and as e grows the oscillations become more 
and more entangled, see Fig. 2(p). 


^ Full synchronization of the subordinates, = 0, and Anti¬ 

phase synchronization of the hub with them, Periodicity 
® All as for F 2 PPI, but chaos instead of periodicity 
® Full Synchronization, Chaos 
Oscillation Death 

Full synchronization of the subordinates, . u = 0, and No 
synchronization of the nub with them, Quasiperiodicity 
Same as FNQ, but chaos instead of quasiperiodicity 



Figure 4. (color online) Longest synchronization windows 
vs. coupling near the boundaries of ranges. (a,b) NSCl, 
£0 = 0.115; (c,d) NSC2, eo = 0.345. and W^_n 

denote the lengths of the longest windows where correspond¬ 
ing d^_„, and vanish. 
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Figure 5. (color online) The fist Lyapunov exponents for 
A = 4, 5, 6 and 7. Data for A = 4 are the same as in Fig. 3(b). 
Vertical lines at e = 0.13 and 0.26 mark an area of wild mul¬ 
tistability where the Lyapunov exponents for different A do 
not coincide. 


C. Stars with A > 4 and wild multistability 

Figure 5 compares first Lyapunov exponents computed 
for stars of different sizes. Observe remarkable coinci¬ 
dence of almost all points, indicating the identity of dy¬ 
namics regardless of N. Additionally we verified it by 
visual inspection of time series and by the analysis of 
synchronization windows. Moreover, below we will pro¬ 
vide rigorous proofs for certain regimes that their ranges 
of stability do not depend on N. 

However, one have also to notice the area where Ai does 
depend on N. This area is marked in Fig. 5 by dashed 
vertical lines. Some of attractors like FPQ (Ai =0) and 
FCP (the lowest branch of Ai) still exists regardless of 
N. But in addition a variety of other attractors appears 
only at certain A, each with a narrow range of stability 
on e. 

Figure 6 illustrates what can be encountered within 
this area. Figure 6(a) demonstrates a periodic regime 
with Ai = —0.136, emerging at A = 5, and e = 0.13. 
Here the subordinates are fully synchronized by pairs: 
d |-5 = 0, and ci |_4 = 0. These pairs in turn are dip- 
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Figure 6. (color online) Examples of dynamics of the star-like 
networks with > 4 within the range of wild multistability, 
(a) TV = 5, e = 0.13 (b) iV = 7, e = 0.24 


natural frequency [14, 26]. 

In our case there are regimes with fully synchronous 
subordinates and the hub not coinciding with them. 
These are FPQ, DPP, FCP, all F 2 PPS, all FAPs, and 
F 2 PC, see Fig. 2. This behavior though reminds remote 
synchronization, can not be classified like this since ac¬ 
tually the hub is also synchronized with the rest of nodes 
with a phase shift. 

The remote synchronization is nevertheless observed 
for our system, see FNQ and FNC in Figs. 2(m) and 
(o,p), respectively. These regimes essentially differ from 
those reported in Ref. [26]. First, the subordinates are 
fully synchronized, while in Ref. [26] phase synchroniza¬ 
tion is reported. Second, in our case oscillators are iden¬ 
tical, but the hub still does not get synchronized with 
the subordinates. Third, FNQ and FNC are quasiperi- 
odic and chaotic regimes, respectively, contrary to the 
periodic case reported in [26] . For these regimes the dif¬ 
ference between the hub and subordinates is found to be 
not necessary to prevent their synchronization. 


dip synchronized with each other, ^^3,4 5 = 0 . The 
hub is in peak-dip synchronization with all subordinates, 
^ 1-2 345 =0- Figure 6(b) represents a chaotic regime at 
N = ’7 and e = 0.24 with Ai = 0.00206. Four of six sub¬ 
ordinates, namely 4,5,6, and 7, are fully synchronized. 
Two others, 2 and 3, are dip-dip synchronized with each 
other. This pair is in turn in peak-dip synchronization 
with the other four. The hub is dip-dip synchronized 
with the nodes 2 and 3 and peak-dip synchronized with 
the other four. 

Altogether, within the area delimited by the dashed 
lines in Fig. 5 there is the large number of attractors with 
highly interwoven basins and narrow stability ranges. As 
a result one can hardly predict the qualitative nature of 
the expected dynamics since each perturbation to e or 
initial condition results in quite different behavior. Sim¬ 
ilar situation was already reported for systems with high 
dimensional phase space, in particular, for coupled map 
lattices, see the review in Ref. [22] . We suggest to re¬ 
fer to it as wild multistability to emphasize its amazing 
richness. 

One more similar area is located in Fig. 5 between e « 
0.8 and e « 0.9. Though the Lyapunov exponents do not 
depend on N here, the basins of coexisting attractors can 
be colocated in complicated manner impeding prediction 
of the dynamics. 


D. Remote synchronization 

As we already mentioned above, remote synchroniza¬ 
tion occurs between two or more nodes without direct 
connections, but linked with the common hub node. The 
important point here is that the hub remains not syn¬ 
chronized with them, that is achieved by detuning its 


V. STABILITY OF CERTAIN REGIMES 
A. Full chaotic synchronization (FSC) 

Stability of the fully synchronized state can be ana¬ 
lyzed using so called Master Stability Function (MSF) [1, 
24, 30] . First we need the Jacobian matrix of the net¬ 
work (1). It has a block form being composed of N x N 
matrices: 


J(t) 


(m A 

V/3I Oj’ 


( 8 ) 


where 

F(t) = -2G(t)[(l-e)I + eK-iA], 

G{t) = diag{a:„ -I- e/i„}, K = diag{A:„}, 

and I is the identity matrix [18] . 

Jacobian matrix at Xn{t) —^ x(t) governs the evolu¬ 
tion of tangent perturbation to the synchronization man¬ 
ifold, where x{t) is produced by a local isolated oscillator. 
Decomposing the perturbations over eigenvectors of the 
matrix (K~^A) one obtains the following equations for 
perturbation amplitudes: 

dp{t + 1) = -v[2x(t)Sp{t)] + Sq{t), 

6q{t +1) = j3 dp, 


where 


i 2 =l + e{(l)-l), (11) 

and (j) is an eigenvalue of (K“^A). We can treat v a free 
parameter and define MSF as a conditional Lyapunov 
exponent, i.e., an average rate of exponential growth of a 
solution of Eq. (10) when x{t) runs along a trajectory of 
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Figure 7. (color online) MSF for the network (1), (2). Vertical 
dotted lines delimit an area where MSF is negative. 


the local oscillator. The graph of MSF is shown in Fig. 7. 
Note that MSF is negative at \v\ < 0.652. 

Using Eq. (2), one can rewrite matrix (K“^A) as 


j^-i^ _ ( ^ 1,1 

\dM,i zm,m , 


and then compute its eigenvalues: 

4>l = 1; = —1, 4'3 = 4'‘l = ■ ■ ■ = (t’N = 0 


( 12 ) 


(13) 


The eigenvector corresponding to (pi has identical ele¬ 
ments being responsible for chaotic dynamics on the syn¬ 
chronization manifold. All other eigenvectors and eigen¬ 
values correspond to transverse perturbations whose van¬ 
ish is the necessary condition of stability of the full syn¬ 
chronization. 

Substituting (f> = —1 and ^ = 0 to Eq. (11), one can 
see that MSF is negative for all transverse perturbations 
when 


e = epsc e [0.348,0.826]. (14) 

Withing this range the full synchronization attractor 
is transversely stable on average. For attractors with 
regular structure this is also a sufficient condition, but 
when the dynamics is chaotic, the synchronization can 
be destabilized by a small noise even when MSF is neg¬ 
ative. This is due to the presence of transversally unsta¬ 
ble invariant sets (cycles, in particular) embedded into 
the synchronization manifold, see [31]. In our case how¬ 
ever the influence of noise (3) is indistinguishably weak. 
One can see in Fig. 3(b) that the first Lyapunov expo¬ 
nent attains the level corresponding to this regime almost 
exactly within the range (14). 

Note that the stability range for the full chaotic syn¬ 
chronization does not depend on N. The reason is that 
the matrix (K^^A) has only three different eigenvalues 
1,-1, and 0 regardless of N. 


B. Oscillations death (OSD) 


The solution of these equations can be expressed via roots 
and ^2 of the polynomial 

{{P — 1)^ — 4a)e^ -|- a(4e — 1) 




(2e-l)4 


= 0. (16) 


There are two couples, Xa = Ci, Xi, = ^ 2 , and Xa = ^ 2 , 
= Cl) corresponding to two forms of OSD. 

Stability of the OSD solution is determined by eigen¬ 
values of the Jacobian matrix (8) at (xa, Xh). The 
eigenvalue problem for is reduced to 

= (17) 


where / is an eigenvalue of F, see Eq. (9). One can see 
that each / produces a couple of /i. Computing / we 
employ the block representation of A (2): 


Y = ({Kae/M)di ,M 

y Kbidnf,! Kfc(l — e)eM,M 

Here Ha = Xa + e(x6 — Xa) corresponds to the hub node 
and Kb = Xb + e{xa — Xb) comes from the subordinates. 
The eigenproblem for F results in 

[/ -I- 2Ka(l - e)][/ -I- 2Kb{l - e)]M = AKaKbt^S. (19) 

Here 5 is one of eigenvalues of djvr.M, 

Ji = M, J 2 = <^3 = ■ ■ ■ = 5m = 0. (20) 

Each J = 0 produces a couple of solutions of Eq. (19). 
However, / = —2 «:q( 1 —e) is spurious that can be verified 
by substituting it to the eigenproblem for F. Thus F 
has the eigenvalue / = —2Kb{l — e) with the multiplicity 
M — 1. Two more eigenvalues correspond to J = M. 
They can be found from Eq. (19) by canceling M. Note 
that none of / depends on the number of the network 
nodes. 

To obtain the stability range for OSD one has to vary e, 
compute /, and And the absolute values of corresponding 
eigenvalues /r from Eq. (17). The regime is stable when 
\^\ < 1. Figure 8 shows the plots of the largest by mag¬ 
nitude \fi\. In the middle area the curves vanish because 
Eq. (16) has complex roots. Thus the range of stability 
computed in this way is 

e = eosD e [0.757,0.856]. (21) 

It agrees well with the range obtained via straightforward 
simulations, see Tab. I. 



C. Full synchronization of the subordinates, and 
complementary synchronization of the hub with 
them (FCP) 


Let Xa and Xf, be states of the hub and the subordinates 
at OSD, respectively. Using Eqs. (1), one can write: 


Xa = a - [Xa + e{xb - Xa)f J- l3Xa, 
Xb = a - [xb + e(xa - Xb)]'^ + Pxb- 


(15) 


In this regime the system oscillates with period two 
so that the hub and the coinciding subordinates switch 
between Xa and Xb- 

Xa = a-[xb + e(Xa - Xb)f + (3Xa, 

Xb = a - [Xa + e{Xb - Xa)Y + I5xb- 


( 22 ) 










Figure 8. (color online) The largest by magnitude eigenvalue 
of the Jacobian matrix (8), (9) for OSD. Vertical lines delimit 
the range of stability. 


Ijll -1-1-1-r 



0 0.143 0.244 0.4 0.6 0.8 e 


Figure 9. (color online) The largest by magnitude eigenvalue 
of the Jacobian matrix for the regime FCP. Vertical lines de¬ 
limit the range of stability. 

Solutions to this equation set can be found as roots of 
the polynomial 

2 - 1 ^ , (l-e)[(l-/3)^(l-e)+4ae]-a ^ 

^ ( 2 e-l )4 

(23) 

Stability of FCP is determined by the eigenvalues /i of 
the matrix (J{,Ja) where Jq is computed as in Sec. VB, 
and Jfe is similar on interchanging Xa and Xb'- 

= (24) 

where / is an eigenvalue of the matrix product (Ft,Fa). 

Using Eq. (18) one can write the matrix (F^Fq) explic¬ 
itly and find that its eigenproblem results in the equation 

(/ ~ ail)(/ ~ b22)M = ( 022 / — 011022 + 0i202i)i5, (25) 

where on = Kb[Ka{l-eY+Hbe\ 012 = Kb{Ka+Kb)e{l-e), 
O2I = KaiKa + Kb)e{l-e), O22 = , &22 = ^^^^(l-e)^, 

and 6 is an eigenvalue of d^.M, see Eq. (20). The root 
/ = On at 5 = 0 is spurious that can be tested by a 
direct substitution. Thus there is an eigenvalue / = 622 
with the multiplicity M —1 and two more eigenvalues are 
the solutions of Eq. (25) at ^ = M. Notice, that as in 
Sec. VB, M is canceled so that the stability range again 
does depend of the network size. 

Varying e we can compute / and then find the largest 
1^1. The plot is shown in Fig. 9. The FCP regime is 
stable where |^| < 1. This range is 

e = epcp e [0.143,0.244]. (26) 


VI. OUTLINE AND CONCLUSIONS 

We considered a variety of dynamics of Henon map 
networks with star-like topology. This is found to be 


amazingly rich. In brief, as the coupling strength grows 
from zero to one the following areas are observed. 

• Non-synchronized oscillations, intermittency of 
synchronization windows. 

• Quasiperiodicity. 

• Wild multistability. Small variations of initial con¬ 
ditions result in periodic, quasiperiodic and chaotic 
solutions. The zoo of regimes is very sensitive both 
to variations of the coupling strength and to the 
network size. 

• Again non-synchronized oscillations, intermittency 
of synchronization windows. 

• Full chaotic synchronization of all network nodes. 

• Coexistence of the full chaotic synchronization, os¬ 
cillation death and periodic oscillations when the 
hub oscillates in anti-phase with fully synchronized 
subordinates. 

• Oscillation death. 

• Remote synchronization of quasiperiodic oscilla¬ 
tions. At narrow range of coupling values it coexists 
with another version of the regime of periodic os¬ 
cillations when the hub is in anti-phase with fully 
synchronized subordinates. 

• Remote synchronization of chaotic oscillations. 

We conjecture that the list of regimes and correspond¬ 
ing stability ranges in most cases remain the same re¬ 
gardless of N. This is checked numerically for V < 7 
and rigorously proved for certain regimes. The exception 
is the wild multistability area where the dynamics is very 
sensitive to variations of N. 

The considered networks demonstrates new examples 
of remote synchronization. These are quasiperiodic and 
chaotic regimes. For these regimes the non-identity of the 
hub and the subordinate nodes is not required to prevent 
their synchronization. 

The most interesting situation occurs in the wild mul¬ 
tistability area. Similar type of dynamics was already 
reported for other systems, in particular for coupled map 
lattices. We suggest to refer to it as wild multistability 
because of its complexity and richness. In fact one can 
hardly predict here even a qualitative nature of the ex¬ 
pected solution. Each deviation in initial conditions and 
coupling strength value is found to result in a new char¬ 
acter of dynamics. This can be either periodic, quasiperi¬ 
odic, or chaotic regimes. Obviously a lot of interesting 
problems arise in this connection. Since star-like motifs 
are the basic building blocks of scale-free networks, the 
detailed study of wild multistability area is very impor¬ 
tant for further understanding of the dynamics of com¬ 
plex scale-free networks. 













9 


ACKNOWLEDGMENTS 

This work was partially supported by a grant of the 
President of the Russian Federation for leading scientific 


schools NSH-1726.2014.2 “Fundamental problems of non¬ 
linear dynamics and their applications” 


[1] S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, and 
D.-U. Hwang, Physics Reports 424 , 175 (2006). 

[2] X. F. Wang, International Journal of Bifurcation and 
Chaos 12 , 885 (2002). 

[3] A.-L. Barabasi, R. Albert, and H. Jeong, Physica A 281 , 
69 (2000). 

[4] A. Arenas, A. Dlaz-Guilera, J. Kurths, Y. Moreno, and 
C. Zhou, Physics Reports 469 , 93 (2008). 

[5] G. V. Osipov, J. Kurths, and C. Zhou, Synchronization 
in oscillatory networks (Springer Science & Business Me¬ 
dia, 2007). 

[6] M. Golubitsky and I. Stewart, Chaos: An Interdisci¬ 
plinary Journal of Nonlinear Science 25 , 097612 (2015). 

[7] I. Belykh, M. Hasler, M. Lauret, and H. Nijmeljer, In¬ 
ternational Journal of Bifurcation and Chaos 15 , 3423 
(2005). 

[8] A. Arenas, A. Dlaz-Guilera, and C. J. Perez-Vicente, 
Phys. Rev. Lett. 96, 114102 (2006). 

[9] X. F. Wang and G. Ghen, Gircuits and Systems I: Fun¬ 
damental Theory and Applications, IEEE Transactions 
on 49 , 54 (2002). 

[10] J. Fan and X. F. Wang, Physica A: Statistical Mechanics 
and its Applications 349, 443 (2005). 

[11] S. Jalan and R. E. Amritkar, Phys. Rev. Lett. 90, 014101 
(2003). 

]12] S. Jalan, R. E. Amritkar, and G.-K. Hu, Phys. Rev. E 
72 , 016211 (2005). 

]13] S. Jalan, A. Singh, S. Acharyya, and J. Kurths, Phys. 
Rev. E 91, 022901 (2015). 

]14] L. V. Gambuzza, A. Gardillo, A. Fiasconaro, L. F. J. 
Gomez-Gardenes, and M. Frasca, Ghaos 23 , 043103 
(2013). 


[15] J. Wang and Y. Zhang, Physics Letters A 374 , 1464 

( 2010 ). 

[16] X. F. Wang and G. Ghen, Physica A: Statistical Mechan¬ 
ics and its Applications 310 , 521 (2002). 

[17] P. V. Kuptsov and U. Parlitz, J. Nonlinear Sci. 22 , 727 

( 2012 ). 

[18] P. V. Kuptsov and A. V. Kuptsova, Phys. Rev. E 90 , 
032901 (2014). 

[19] A. Roxin, N. Brunei, and D. Hansel, Phys. Rev. Lett. 
94 , 238103 (2005). 

[20] D. Angeli, J. E. Ferrell, and E. D. Sontag, Proceedings 
of the National Academy of Sciences of the United States 
of America 101 , 1822 (2004). 

[21] D. Angeli, European Journal of Control 15, 398 (2009). 

[22] U. Feudel, International Journal of Bifurcation and 
Chaos 18 , 1607 (2008). 

[23] A. N. Pisarchik and U. Feudel, Physics Reports 540, 167 
(2014), control of multistability. 

[24] L. M. Pecora, Phys. Rev. E 58 , 347 (1998). 

[25] Z. Ma, G. Zhang, Y. Wang, and Z. Liu, Journal of 
Physics A: Mathematical and Theoretical 41, 155101 
(2008). 

[26] A. Bergner, M. Frasca, G. Sciuto, A. Buscarino, E. J. 
Ngamga, L. Fortuna, and J. Kurths, Phys. Rev. E 85, 
026208 (2012). 

[27] M. Henon, Communications in Mathematical Physics 50, 
69 (1976). 

[28] A. Politi and A. Torcini, Chaos 2 , 293 (1992). 

[29] U. Feudel, S. Kuznetsov, and A. Pikovsky, Strange non- 
chaotic attractors (World Scientific, 2006). 

[30] L. M. Pecora and T. L. Carroll, Phys. Rev. Lett. 80, 2109 
(1998). 

[31] J. Milnor, in The Theory of Chaotic Attractors (Springer, 
2004) pp. 243-264. 



